Axisymmetric versus Non-axisymmetric Vortices in Spinor Bose-Einstein Condensates 
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The structure and stability of various vortices in F — 1 spinor Bose-Einstein condensates are 
investigated by solving the extended Gross-Pitaevskii equation under rotation. We perform an ex- 
tensive search for stable vortices, considering both axisymmetric and non-axisymmetric vortices and 
covering a wide range of ferromagnetic and antiferromagnetic interactions. The topological defect 
called Mermin-Ho (Anderson- Toulouse) vortex is shown to be stable for ferromagnetic case. The 
phase diagram is established in a plane of external rotation fi vs total magnetization M by compar- 
ing the free energies of possible vortices. It is shown that there are qualitative differences between 
axisymmetric and non-axisymmetric vortices which are manifested in the fi- and M-dependences. 

PACS numbers: 03.75.Fi, 67.57.Fg, 05.30.Jp 
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I. INTRODUCTION 

The experimental achievement of Bose-Einstein con- 
densation (BEC) in the trapped atomic clouds ||, |^| 
has opened up a novel field to investigate fundamen- 
tal problems such as the relationship between the su- 
perfluidity and BEC. Owing to recent advances of ex- 
perimental techniques, several groups have succeeded in 
creating quantized vortices with various procedures in 
the magnetically trapped BEC{|, §, || % §, §], where 
the condensate is described by a scalar order parameter. 
Furthermore, the atomic gases with the hyperfine spin 
F = 1 called "spinor BEC" have been Bose-condensed 
via optical methods [^0[ [llj which can keep the atomic 
"spin" states degenerate and active. As shown recently 
by Klausen et al, [jl2|, the spin-dependent interaction of 
two 87 Rb atoms is ferromagnetic. Thus, we now have 
spinor BEC's with both antiferromagnetic ( 23 Na)[jllJ and 
ferromagnetic interactions. 

Such scalar and spinor BEC's are analogous to super- 
fluid 3 He and 4 He. These superfluid Heliums, however, 
have rather strong interactions. Indeed, the conden- 
sate fraction in superfluid 4 He is only 10% of the total. 
By contrast, BEC of the atomic gases have advantages 
for both theoretical and experimental treatments due to 
their weak interactions; here almost all the atoms are 
able to be Bose-condensed. It is possible to directly ob- 
serve dynamical behaviors of the condensate with optical 
methods, providing us an opportunity to quantitatively 
investigate the new quantum fluid p^]. 

The standard Hamiltonian for the spinor BEC have 
been introduced by Ohmi and Machida [Q, and Ho[p^[, 
who pointed out the richness of the exotic topological de- 
fects. Topological structures, such as skyrmion, meron, 
Mermin-Ho (Anderson- Toulouse) texture, and monopole, 
play an important role in various fields of physics. They 
provide a common framework to connect diverse field, 



thereby enhancing mutual understanding^, p7{ . Spe- 
cially, since F = 1 ferromagnetic BEC can be described 
by order parameters similar to the superfluid 3 He-A, the 
coreless Mermin-Ho vortex may be favored in the ferro- 
magnetic BEC@|l9 . 
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Similar topological structures, called skyrmion in gen- 
eral, have been proposed in the spinor BEC. Al Khawaja 
and Stoof [^0| studied a skyrmion in the F = 1 ferromag- 
netic BEC and concluded that it is not a thermodynam- 
ically stable object without rotation. By considering the 
effect of the external rotation, however, we have shown 
recently|2l| that this topological defect can be stable. In- 
genious proposals have been madcj22|, 53, [2^] on how 
to create it and detect it. Yip[^6| has performed a sys- 
tematic study on vortex structures and presented several 
axisymmetric and non-axisymmetric vortices for F = 1 
antiferromagnetic BEC. Recently, Isoshima et al. p7j, |2^j 
have carried out an extensive study of axisymmetric vor- 
tices to provide a vortex phase diagram in a plane of the 
rotation and the magnetization for both the antiferro- 
magnetic and ferromagnetic cases. 

In this paper we examine the stability of various vor- 
tices for both the F — 1 antiferromagnetic and ferro- 
magnetic BEC trapped in a two-dimensional harmonic 
potential. We have removed the previously imposed re- 
striction in the axisymmetric case that winding numbers 
are less than or equal to unity. The continuous vortices 
such as the Mermin-Ho vortex will also be shown to be 
favored over the singular onesp7|, p8| and the other non- 
axisymmetric ones [Em. We demonstrate the stability of 
such vortices and discuss differences between the axisym- 
metric and non-axisymmetric configurations. We also de- 
termine the vortex phase diagrams for the antiferromag- 
netic and ferromagnetic cases. By comparing the relative 
free energies of the possible vortex configurations and the 
phase-separated state which may occur in the ferromag- 
netic situation, the validity of assuming uniformity along 
z direction will be checked. 

This paper is organized as follows. In Sec. II, we first 
present the extended Gross-Pitaevskii equation for the 
spinor BEC, and then explain the numerical procedure 
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to find local minima of the energy functional. Section III 
enumerates possible vortices for axisymmetric and non- 
axisymmetric cases, and then discuss the stability of the 
Mcrmin-Ho vortex. Section IV presents the phase dia- 
gram for the ferromagnetic and antiferromagnetic cases 
in the plane of external rotation Q vs total magnetiza- 
tion M obtained by comparing the free energies. Here, 
the qualitative differences between the axisymmetric and 
non-axisymmetric vortices are discussed by showing the 
il- and M-dependence. The final section is devoted to 
conclusions and discussions. 



II. THEORETICAL FORMULATION 

A. Extended Gross-Pitaevskii equation 

We consider Bose condensed F=l spinor BEC's with 
internal degrees of freedom for both ferromagnetic and 
antiferromagnetic cases. Here the order parameters are 
characterized by the hyperfine sublevels mp = 1,0,-1. 
We start with the standard Hamiltonian by Ohmi and 
MachidalHl, and HolTl: 



7i, 



H - J drfl ■ *}(r x P)*i 



n 



dr 



a ijkl 



Here 



h(v) 



n 2 v 2 

2m 



V(v) 



(1) 



(2) 



is one-body Hamiltonian. The quantity V(r) — 
im(27rtv) 2 {x 2 + y 2 ) is the external confinement poten- 
tial such as an optical potential. The scattering lengths 
ao and a-i characterize collisions between atoms through 
the total spin and 2 channels, respectively, g n = 
At^_ . a +2a 2 j g interaction strength through the "den- 
sity" channel, and g s = ■ is that through 
the "spin" channel. The subscripts a = (x,y,z) and 
i,j,k,l = (0, ±1) correspond to the above three species. 
The chemical potentials for the three components /Ltj 
(i = 0,±1) satisfy fix — /io = //o — M-i- We introduce 
fi = no and fj,' = \i\ — /Jo- The angular momentum oper- 
ators F a (a = x, y, z) can be expressed in matrices as 



F, 




0-10 
Fy = -}= I 1 -1 
v2 \ o 1 




(3) 



Following the standard procedure, the extended Gross- 
Pitaevskii (GP) equation in rotation frame is obtained as 



[{h(r)-m + gnJ2\i>i\ 2 }Sij 



a kl 

-ihU ■ V x rSijtyj = 0. (4) 

These coupled equations for the j-th condensate wave 
function ipj = (tyj) (j — 0, ±1) are used to calculate var- 
ious properties of vortices in the following. Here we take 
the external rotation as f2 = £lz and assume uniformity 
along z direction. 



B. Numerical procedure 

The stationary states of the extended GP equation are 
defined as local minima of the energy functional 



where and E s are defined by 
9. 



fi'M-Sl • L,(5) 



a I fej 



f2 • L denotes 



ft L = -im J drY^ip* 



X—- —)r 

dy dx } 



(6) 
(7) 

(8) 



and M — J dv^2 i i\tfj i (r)\ 2 is the total magnetization. 

The numerical algorithm used to minimize the energy 
functional in scalar BEC|29j ||l], |3^] can be extended 
to the present system. Following this procedure, the ini- 
tial tpj given randomly are modified using the local gra- 
dient of the energy functional E[ijji,ip*] as 



i>3 ( T + Ar ) = W r ) ^ At. 



(9) 



This equation means that the order parameters tpj, pa- 
rameterized by a 'fictitious time' r, roll along the slope 



3 



of the energy functional. Equation (^) is rewritten as 



<h(r)-(ji(T)+»'j)+g n J2\i>k 

. { k 



Ojk 



dy 



y 



Ox 



Ojk 



(10) 



This is the Gross-Pitaevskii equation for imaginary times 
t = it. In each time step, /Lt(r) is adjusted to preserve 
the total number of particles in the system 



N 



E / d ^M' 



(ii) 



For t — ► oo , tpj converges to the stationary state, corre- 
sponding to one of the local minima of the energy func- 
tional (||). For t = oo tpj satisfies 



= o, 



(12) 



and Eq. (n0|) becomes equivalent to the extended GP 
equation (0). 

We take the initial state of each component as 



i) j (r,T = 0) = \Ai T F(r) rjj exp[i5j(r)], 



(13) 



where titf^) is the density profile within the Thomas- 
Fermi (TF) approximation: 



TlTF(r) = 



i^^f^- for g s > 

tZJLzY. for g s < ' 



(14) 



and rjj represents the ratio of each component. The phase 
is given by 

.,( fc )fl( fc ) _■_ 



E w i 



w 



(fc) 



(15) 



is the winding number of the j'-th condensate, 6^ is 
the polar angle of the coordinate (x^ , j/W ) whose origin 
is located at the fc-th vortex core, and aj is relative phase 
between the three components. 

It is convenient to describe the condensates in terms 
of the three components ip a (a = x, y, z) where the quan- 
tization axis is taken along the a direction: 



i /-I 1 
^(r) = _ -i -i 
^(r) I V2 \ V2 



Mr) | (16) 
lM(r). 



We then define a couple of real vectors as 

m= (m x ,m y ,m z ) = Re(ipa;,ipy,ip z ), (17) 
n = (n x ,n v ,n z ) =1m.(ip x ,ip y ,ip z ). (18) 

The Z-vector, which points the direction of the local mag- 
netization, is defined as I = m x n. The corresponding 
unit vector is denoted by I. 



C. Calculated system 

The actual calculations are carried out by discretizing 
the two-dimensional space into 51x51 mesh. We have 
performed extensive search to find stable vortices, start- 
ing with various vortex configurations, covering a wide 
range of the ferromagnetic and the antiferromagnetic in- 
teraction strength, g s /g n = —0.2 ~ 0.2, and examin- 
ing various axisymmetric and non-axisymmetric vortices. 
See Ref . |27], H for the classification of possible vortices 
in the axisymmetric case. We use the following param- 
eters: the mass of a 87 Rb atom m=1.44x 10 _25 kg, the 
trapping frequency ^ r =200Hz, and the particle number 
per unit length along the z axis n z =2.0x 10 3 /^m. The 
results displayed here are for g s /g n — —0.02 (ferromag- 
netic case) and g s /g n — 0.02 (antiferromagnetic case). 
The external rotation frequency Q is normalized by the 
harmonic trap frequency. 



III. VORTEX STRUCTURE 

The vortex configurations are characterized by the 
combination of the winding number Wj of tpj (j — 0, ±1) 
denoted by (wi, Wq, W-i), where Wj denotes the phase 
change by litWj when the wave function goes around the 
phase singularity. 

The spin term (]?]) of the total energy is rewritten as 



n»- |2^(rty_i(r)-^(r)| 



(19) 



where n(r) = l^'MI 2 is the total density. The spin 
texture in the ground state without rotation is deter- 
mined by this energy. By minimizing Eq.([l9|), the rela- 
tive phases aj are shown to satisfy 



2uq = c±\ + a_i 



(20) 



corrcspon 



ids 



where n is an integer, and the odd (even) n 
to the antiferromagnetic (ferromagnetic) situation 



A. Axisymmetric vortex 

It also follows from Eq. ( |l9|) that (wi, wq, W-i) of ax- 
isymmetric vortices satisfies 



2wq = w\ + w- 



(21) 



Thus, possible candidates for the stable state are the 
non-vortex state (0, 0, 0) and the vortex configurations: 
(1,0, —1), (1, x, 0), (0, x, 1), and (1, 1, 1), which exhaust 
all the combinations of the winding numbers less than or 
equal to unity. The thermodynamic stability of these vor- 
tices are demonstrated in Ref.|2^]. Here we concentrate 
on the possibility of combinations with higher winding 
numbers, which will be shown to be stable only in the 
ferromagnetic case. 
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Figure 1 displays the density and /-vector profiles of a 
new continuous vortex (0,1,2) found stable for the fer- 
romagnetic interaction (g s < 0) . It is seen that tpi with 
zero winding number w\ — occupies the central region 
of the harmonic trap and ■0-1 with the higher winding 
number w-\ = 2 fill in the circumference region. The 
intermediate region is occupied by tpo component which 
has a singularity wq = 1 at the center of the trap. The 
resulting total density is non-singular and have a smooth 
spatial variation described by a Gaussian form except 
for the outermost region. This vortex is equivalent to 
the topological structure called the Mermin-Ho vortex in 
superfluid 3 He|ll| and Skyrmion in general ]2C||. 

The axisymmetric (0, 1, 2) vortex may be parametrized 

as 



-01 (r) 
V'o(r) 
0-i(r) 



cos 



2 £ 

I, 



= \J n(r) | v / 2e 4v sin ^ cos ^ 
e 2 ^ sin 2 § 



(22) 



where the bending angle (3{r) runs over < (3{r) < n and 
(j> signifies the polar angle in polar coordinates. The spin 
direction is denoted by the /-vector and is given as l(r) = 
zcos/3 + sin/3(cos0x + sin0y) where (3 varies from {3(0) = 
to (3(R) = § (= tt) for MH (Anderson-Toulouse (AT)) 
(R is the outer boundary of the cloud). Thus the spin 
moment is flared out to the radial direction and at the 




-3.0 -2.0 



-1.0 1.0 2.0 3.0 4.0 
x(y=0) [/im] 




circumference it points outward for MH and downwards 
for AT (for schematic /-vector structure, see Fig. 18 in 
Ref.@ ). 

In Fig.|| we show the spatial dependence of the l z - 
component along the radial direction, namely, the spa- 
tial dependence of the bending angle f3(r) for the MH 
vortex. As the magnetization M decreases, the local 
magnetization in the condensate surface changes from 
positive to negative passing through zero. It means that 
the /-vector in this vortex flares out radially to orient al- 
most horizontally (3(r = R) = \ for M/N ~ 0.5 and to 
point downward for /3(r = R) = tt for M/N ~ 0. The 
former (latter) corresponds literally to the Mermin-Ho 
(Anderson- Toulouse) vortex. This is simply because as 
M decreases, the spin-down component 0_i with w — 2 
increases in the outer region. Thus we can control these 
MH and AT vortices by merely changing the total mag- 
netization. 

As pointed out in the previous paper pl[| , however, the 
situation is completely different from the case of super- 
fluid 3 He- A where the stability of the MH vortex is due to 
the constraint that the /-vector be perpendicular to the 
vessel wailful. These vortex configurations in ferromag- 
netic BEC are created naturally under the condition of 
a given total number and magnetization, both of which 
are well controlled in a harmonic trap experiment. 

In comparison with the (0, 1,2) vortex, the spin tex- 
tures of other vortex configurations, such as the (1, x , 0) 
and (1,0,-1) vortices, have a different nature fl27f . In 
(l,x,0) vortex, the spin moment is suddenly reversed 
near the vortex core of ipi component because of the ab- 
sence of 0o component. In (1,0,-1) vortex it can vary 
continuously around the vortex core. In this configura- 
tion, however, since the condensate at the center of the 
trap consists only of the polar state the spin texture 
has a singularity. Thus only (0, 1, 2) vortex can have a 
non-singular and continuous spin texture under slow ro- 
tation. 

It is easy to calculate the total angular momentum L z 
of the axisymmetric vortices; by using the total particle 
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FIG. 1: Properties of the (0, 1,2) vortex at D. = 0.35 and 
M/iV = 0.21. (a) Density Profile; (b) and l y ; (c) density 
map of l z . The bold line in (a) denotes the total density n(r) 
and the thin lines show the density of each component \ipj\ 2 - 



FIG. 2: Spatial dependence of the Z z -component along the 
radial direction at M/N=0, 0.46, and tt = 0.37. The dashed 
line shows cos/9(r) with the bending angle f3(r) = nr/R (R — 
2.85fim) for the AT vortex. 
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.4 -3 -2 -1 1 2 3 4 -4 -3 -2 -1 1 2 3 4 
x(y=0) [jum] 

FIG. 3: Density profiles of the (0, 1,2} vortex in an antifer- 
romagnetic interaction of g a /g n = 0.005 (on the left) and the 
non-magnetic interaction of gs/gn = (on the right) situa- 
tion. 



number and the total magnetization, it is simply written 

as 



hN 



wa+w —, 



(23) 



where the total magnetization is written as M = 
J dvn{v)l z and we have introduced Wj — wq — jw' . Thus 
the spin textures with a net spin polarization carry the 
angular momentum, i.e. the superflow. For (0,1,2) vor- 
tex, J*. = 1 — M> . This simple formula has the following 
physical meaning. (i)At M = N, L z — because ipi has 
no winding. (ii)At 4$ = |, is exactly equal to Ti/2, 
corresponding to the MH vortex. 

In the (0, 1, 2) state with the higher winding, the non- 
winding component ipi works as a "pinning potential" 
for the remaining tpo and ip—i, thereby making the state 
stable in the lower rotation drive. In particular, with 
w-i = 2 is stabilized by the presence of the ipi due to the 
ferromagnetic interaction. For a very small antiferromag- 
netic interaction (g s = 0.005c/ n ) and non-magnetic case 
(Ss/ffn = 0), the vortex with the W-% — 2 becomes unsta- 
ble and splits into a couple of id_i = 1 vortices (see Fig.^J). 
This configuration is equivalent to the vortex found by 
Yip (phase IV in Ref. ) an d is always unstable for the 
large g s /g n (> 0). 



B. Non-Axisymmetric vortex 

To investigate the possibility of non-axisymmetric vor- 
tices, let us first recapitulate the axisymmetric (1,1,1) 
vortex. In the axisymmetric case, the total density of 
the (1, 1, 1) vortex is equivalent to the one in scalar BEC, 
having the singularity at the center of the trap where the 
potential energy is minimum. The axisymmetric singu- 
lar (1, 1,1) vortex is always unstable even in the higher 
rotationp8|. However, by displacing the vortex cores of 
each component from the center of the trap, the (1, 1, 1) 
vortex can be stable as a non-axisymmetric non-singular 
type in the lower rotation frequency. This works favor- 
ably to gain more condensation energy at the center of 
the trap. A similar situation is seen near the Q C 2 in the 
homogeneous system, where the singular vortex lattice, 
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(b) 

FIG. 4: Density profiles in the (1, 1, 1} split-(I) (on the left) 
and split-(II) vortex (on the right) at Q — 0.35 and M/N = 
0: (a) Non-magnetic interaction case (<? s /<?n = 0.0) and (b) 
g 3 = 0Ag n . The bold line show the total density n(r) and 
the thin and dashed lines present the density profiles of the 
internal structures, respectively. 



called the Abrikosov lattices, is never favored in the en- 
tire antiferromagnetic region [Q. 

For the antiferromagnetic case, ipi and overlap 
to minimize the spin-dependent energy E s (r). In FigJ^, 
we show the density profiles and <? s -dependence of two 
different (1, 1, 1) vortices. In the non-magnetic limit 
(fl , s/5n = 0), the two states are completely equivalent. As 
the spin interaction g s /g n increases, striking differences 
grow between the two states. For g s /g n ^0A, the vortex 
cores of the state presented in the left of Fig.|| (b) (the 
split- (I) state) collapse and the amplitude of each order 
parameter cannot recover near the vortex cores, i.e. this 
state behaves like the phase separation in x-y plane. In 
contrast, the state shown in the right of Fig.^ (b) (the 
split- (II) state) forms the regular cores. As a result, the 
split- (II) state is energetically favorable over the split- (I) 
state for the antiferromagnetic situation. In very small 
spin interaction range (g s ~ 0.02g„), which is a realis- 
tic parameter, however, the split-(I) state can be favored 
over the split- (II), though the energies of two states are 
very close to each other. 

A third vortex configuration is displayed in FigJ^ where 
vortex cores of each component are displaced to form a 
triangle. In the non-magnetic limit, the three compo- 
nents are completely equivalent and the three singular- 
ities form a regular triangle. As the spin interaction g s 
increases, the ipi and tp-i components overlap locally be- 
cause of the antiferromagnetic interaction and this vor- 
tex configuration starts to deform from a regular triangle 
shape. For larger spin interactions (g s > 0.4), this vortex 
becomes unstable due to the overlap between ipi and ip-i 
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(a) 



(b) 




(c) 

FIG. 5: Density profiles of each component |?/ii(r)| 2 , ?/>o(r)| 2 , 
and |-i/)_i(r)| 2 in (1,1,1) triangle vortex at Q = 0.35 and 
M/N = 0: (a) Non-magnetic interaction case (gs/g-n = 0.0), 
(b) g s = 0.02g n , and (c) g s = 0.2g n . The total density pro- 
file almost agrees with that in the vortex-free state, i.e. the 
non-singular type. 



components. 

Figure |^ (a) shows the density profiles of a non- 
axisymmetric non-singular (1,1,1) vortex for the ferro- 
magnetic case. Two singularities of tpi and t^-i are dis- 
placed from the center. The ipo component with the sin- 
gularity at the center of the trap plays the role to prevent 
the phase separation favored in the ferromagnetic spin in- 
teraction. The spin texture in this state is displayed in 
Figs.p (b) and (c) where the spin moments flip at the 
center of the trap. 

In comparison with axisymmetric types, these non- 
axisymmetric vortices have the advantage that they can 
easily adjust themselves for a change in f2. As fl in- 
creases, the two or three separate singularities adjust 
their mutual distance from the center and change the 
value of L z to gain the energy from the term —QL Z . 
In this sense, the non-axisymmetric vortices are flexible 
against a change in Q compared with the axisymmetric 
ones. 



IV. PHASE DIAGRAM 

The phase diagrams are shown in a plane of the exter- 
nal rotation and the total magnetization by comparing 
the energies of various vortex configurations: 
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FIG. 6: (a)Density profile of the condensates, (b)the l x and 
l y and (c)the density map of l z for the non-axis-symmetric 
(1, 1, 1) vortex inO = 0.35 and M/N=0.0. The bold line is 
the total density n(r) and the thin lines show the density of 
each component \ipj\ 2 ■ 



It is noted that the region fl < 0.38 considered here cor- 
responds to the single vortex region in the scalar BEC 
case|55|. Thus, the present single- vortex consideration 
may also be justified for £1 < 0.38. 

The resulting phase diagram is displayed in FigjT] for 
(a) the ferromagnetic case and (b) the antifcrromagnetic 
case. For the ferromagnetic large area of the 

il — M plane is occupied by the (0, 1, 2) vortex, includ- 
ing MH and AT. The non-axisymmetric (1,1,1) vortex 
and the (1,0,-1) vortex are stabilized near M = and 
M = N, respectively. We find a large empty area in the 
intermediate M/N region where neither single- vortex nor 
vortex-free states are stabilized at all because the phase 
separation in the ferromagnetic case prevents forming a 
uniform mixture of the three components even when the 
circulation is absent in the vortex-free state. 

In the phase diagram for the antiferromagnetic case, in 
contrast, everywhere is occupied by a stable phase. This 
result is consistent with Fig. 2(a) of Ref. |28[ over a wide 
range, except for the presence of two non-axisymmetric 
types near M ~ 0, i.e. (1, 1, 1) split-(I) and triangle vor- 
tex. The (1, 1, 1) triangle vortex is energetically indistin- 
guishable from the phase-I vortex obtained by Yip|26|. 
The phase- IV vortex given in Fig, 3 of Ref. |^6| is found 
unstable for g s = 0.02g„ and does not appear in this phase 
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FIG. 7: Phase diagram for (a)the ferromagnetic state (g s = — 0.02<7„) and (b)the antiferromagnetic case (g a = 0.02q T1 ). The 
dashed line in (a) denotes the boundary where the lowest quasiparticle energy of the (0, 1, 2) vortex becomes negative [[2 lfl. 



diagram. 

It is found for both cases that the stabilities of the non- 
axisymmetric vortices are restricted in a narrow region. 
On the other hand, the axisymmetric vortices are stable 
in a large area. As discussed in Section III, since finite 
density of a component at the vortex cores of the others 
supports its stability, the non-axisymmetric vortex be- 




(a) 



0.2 0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 0.38 
Q. [trap unit] 




0.2 0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 0.38 
Q. [trap unit] 



FIG. 8: ^-dependence of the angular momentum L z /N at 
M/N = 0: (a) the ferromagnetic case (g s = — 0.02g n ) and (b) 
the antiferromagnetic case (g s = 0.02g„). 



comes unstable in increasing M where the number of ipi 
grows while the others shrink. 

Figure || shows the ^-dependence of the angular mo- 
mentum at M = for the two cases. For axisymmetric 
types, as shown in Eq. (f23|) , the angular momentum of 
the system L z is fixed for a given M. Thus there is the 
need of changing the winding combinations (u>i , Wq, W-\) 
so as to increase L z , which means that the axisymmetric 
vortices do not have the adaptability for changes in fl. 



V. DISCUSSION 

The phase-separated state with w — is expected to 
be stable in a large empty region of Fig. 7 (a). In this 
state, and ^-1 components phase-separate along z- 
direction due to the ferromagnetic interaction. Namely, 
an arbitrary x-y cross-section consists of only ipi or 
component, and the spin-polarized state with M/N = ±l 
is piled up along the z-direction. Neglecting the con- 
tribution from the boundary layer associated with the 
phase separation, we can estimate that the energy of this 
phase-separated state with w = and w = 1 is equal to 
the energy of the scalar BEC with w — and w = 1 , 
respectively. 

We compare in Fig.^ (a) the free energies of the (0, 1, 2) 
vortex state and the phase-separated state. As shown in 
Fig.^ (a), the (0, 1, 2) vortex with the three components 
is energetically favored over the phase-separated state, 
where the energy of the phase-separated state is given by 
the energy of the (0, 1,2) vortex at M/N = 1, i.e. the 
scalar BEC with w = 0. Thus, the composite state of the 
three components, such as the MH vortex, becomes "lo- 
cally" stable under a rotation drive while the composite 
state may have "global" stability. It is possible to per- 
form a similar discussion for the phase-separated state in 
higher rotation. The result of Fig.[) (b) shows that the 
(1, 0, —1) vortex is favored over the phase-separated state 
with the winding, whose total density corresponds to the 
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VI. CONCLUSION 

We have presented possible vortex structures and the 
vortex phase diagram in the plane of external rotation 
f2 and the total magnetization M/N for the both cases 
of antiferromagnetic (g s — 0.02<?„) and ferromagnetic 
(g s = — 0.02g n ) interaction. We have investigated the 
thermodynamic stability of the possible vortex configu- 
rations by solving the extended Gross-Pitaevskii equation 
for the spinor BEC with F = 1. 

For the ferromagnetic case, the stability of the con- 
tinuous vortex, called the Mermin-Ho and Anderson- 
Toulouse vortex, is demonstrated, but these topological 
structures are found never stable under no rotation drive. 
Furthermore, these vortices can exist in the intermediate 
process (see Fig. 5 in the Ref . p7| ) proposed by Isoshima 
et al. J|6|, [37j, i.e. it may be created by making use of spin 
texture. 

We have also discovered a couple of new non- 
axisymmetric vortices besides the two vortices found by 
Yip|{26| for the antiferromagnetic case. The conventional 
singular vortex is found to be never favored in spinor 
BEC for both cases. It means that the total density pro- 
file is always non-singular and continuous. Therefore, the 
experimental procedure to image the magnetic patterns 
for each vortex configuration will be required as a special 
technique to identify these vortices. 



FIG. 9: The M-dependence of the energy for the ferromag- 
netic case: (a)the (0,1,2) and (b)the (1,0,-1) vortices. The 
energies of the phase-separated states with w = and w — 1 
correspond to the energies of the (0, 1, 2) vortex at M/N — 1 
and (1,0,-1) vortex at M/N = l, respectively. 



conventional singular vortex in the scalar BEC. 
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